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ABSTRACT 

We investigate the high-redshift quasar luminosity function (QLF) down to an apparent magnitude 
of Iab = 25 in the Cosmic Evolution Survey (COSMOS). Careful analysis of the extensive COSMOS 
photometry and imaging data allows us to identify and remove stellar and low-redshift contaminants, 
enabling a selection that is nearly complete for type-1 quasars at the redshifts of interest. We find 155 
likely quasars at z > 3.1, 39 of which have prior spectroscopic confirmation. We present our sample in 
detail and use these confirmed and likely quasars to compute the rest-frame UV QLF in the redshift 
bins 3.1 < z < 3.5 and 3.5 < z < 5. The space density of faint quasars decreases by roughly a factor 
of four from z ~ 3.2 to z ~ 4, with faint-end slopes of /3 ~ —1.7 at both redshifts. The decline in space 
density of faint optical quasars at z > 3 is similar to what has been found for more luminous optical 
and X-ray quasars. We compare the rest-frame UV luminosity functions found here with the X-ray 
luminosity function at z > 3, and find that they evolve similarly between z ~ 3.2 and z ~ 4; however, 
the different normalizations imply that roughly 75% of X-ray bright active galactic nuclei (AGN) 
at z ~ 3 — 4 are optically obscured. This fraction is higher than found at lower redshift and may 
imply that the obscured, type-2 fraction continues to increase with redshift at least to z ~ 4. Finally, 
the implications of the results derived here for the contribution of quasars to cosmic reionization are 
discussed. 



1. INTRODUCTION 

The evolution of the quasar luminosity function 
with redshift is a key observational constraint on the 
growth of supermassive black holes (SMBHs) over cos- 
mic t i me dRichstone et alJll998t [Kauffmann fc Haehneltl 
l2000t IWvithe fc Loebl 120031: IMarconi et all 120041 ). The 
behavior of the QLF places constraints on the duty 
cycles of quasars, the growth history of SMBHs, and 
the coevolution of black holes and the ir host galax- 
ies ([Hopkins et al.ll2006l lUeda et aL1l2003l and references 
therein) . 

The QLF also determines the cumulative ionizing back- 
ground radiation due to quasars. The faint end of the 
QLF at high redshift is of particular interest in this re- 
gard because faint quasars contribute substantially to 
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the total ionizing background due to quasars. While the 
peak in quasar activity around z ~ 2 — 3 is responsi- 
ble for Hell reioniza tion at z ~ 3 (Reimcrs et alJH997t 
ISokasiane t al. 2002}), the relative contribution of quasars 
to hydrogen reionization at z ~ 6 — 10 is not as well con- 
strained. The observed decline in space density of highly 
luminous quasars at z > 3 has been taken as evidence 
that quasars contribute negligibly to hy drogen reioniza- 
tion (jMadau et al.lll999HFan et al.ll2006l) . However, only 
the most luminous quasars at high redshift can be found 
in surveys such as SDSS, leaving the contribution of faint 
quasars to the ionizing background unknown. 

The faint end is challenging to study at high redshift 
because of the need for survey fields that are both deep 
enough to detect faint sources and wide enough to find a 
statistically significant number of quasars. In addition, 
follow-up spectroscopy is difficult due to the faintness of 
the sources. Two groups hav e recently investigat ed the 
faint end of the QLF at z ~ 4: IGlikman et al.ll20lil (here- 
after Gil) using parts of the Deep Lens Survey (DLS, 



alter uiij using par ts ot tne ueep Dens Survey (UL,t> 
IWittman et aL|l2002h and NO AO D eep Wide-Field Sur 



vey fNDWFS. Uannuzi fc Devlll999D fi elds, with a com- 
bined sky coverage of 3.76 deg 2 , and llkeda et"all l20Tll 
(hereafter 1 11) using the HST- ACS region of the COS- 
MOS field (j Scoville et al. 1 120071) with a sky coverage of 
1.64 deg 2 . The value of the faint-end slope /3 measured in 
Gil, —1.6io'g, is consistent with the value of — 1.671q 17 
reported in 111 for the COSMOS field. 

While these studies agree on the faint-end slope of the 
QLF at z ~ 4, they disagree on the absolute space den- 
sity of low-luminosity quasars by roughly a factor of four, 
with a higher space density reported in Gil. This dis- 
crepancy cannot be attributed to cosmic variance (see 
§ 8.1), and leads to different pictures of quasar evolu- 
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tion. The result reported in Gil implies that the decline 
in the space density of faint quasars with redshift after 
the peak at z ~ 1 — 2 eventually stops and possibly re- 
verses, which could make the contribution of quasars to 
cosmic reionization significant. 

Both 111 and Gil use broad-band optical color selec- 
tion to identify quasars. A significant uncertainty in 
this approach is the selection function, or the fraction 
of quasars that are selected as a function of redshift and 
magnitude. The selection function is usually determined 
with Monte Carlo simulations, which are sensitive to the 
assumed distribution of quasar spectral energy distribu- 
tions (SEDs) in the survey field. This distribution is 
uncertain, and mismatches between the assumed and ac- 
tual distributions can give rise to significant errors in 
completeness estimation and thus the derived luminosity 
function. 

In addition to selection function uncertainty, contam- 
ination from stars and galaxies can be substantial and 
difficult to quantify. Spectroscopy allows a robust de- 
termination of the contamination, but due to the large 
amount of telescope time required, often only a subset 
of the candidates can be observed. Extrapolating to the 
complete sample can result in large errors in the esti- 
mated contaminant fraction. 

Here we take a different approach (described in § 2), 
using the COSMOS survey data to select and identify 
high-redshift quasars with high completeness and low 
contamination. Our principal goals are to (1) determine 
the faint end of the QLF from 3.1 < z < 5 in COSMOS, 
thereby providing an independent check of the result re- 
ported by 111, as well as a direct determination of the 
evolution of the faint-end over this redshift range, (2) 
compare the rest-frame UV quasar luminosity function 
with the rest-frame 2-10 keV X-ray luminosity function 
for sources at z > 3 in order to investigate the obscured 
fraction of bright AGN at these redshifts, and (3) use our 
results to make inferences regarding the evolution of the 
QLF and the likely contribution of quasars to reioniza- 
tion, both of Hell at z ~ 3 and HI at z ~ 6 — 10. 

The structure of this paper is as follows. In § 2 we 
provide an overview of the method used to determine 
the QLF. In § 3 we present our initial quasar selection. 
In § 4 we describe the tests used to verify the high com- 
pleteness of the selection. In § 5 we describe how we 
separate stellar and low-z contaminants from the high- 
redshift quasar population of interest. In § 6 we present 
the final sample of confirmed and likely quasars and dis- 
cuss the X-ray properties of the sample. In § 7 we use 
the final list of likely quasars to compute the luminos- 
ity function at z ~ 3.2 and z ~ 4. In § 8 we discuss 
sources of error that may influence our results. In § 9 we 
compare our results with the X-ray luminosity function 
at similar redshifts and discuss the implications for the 
obscured quasar fraction at these redshifts. In § 10 we 
discuss the ionizing background due to quasars. In § 11 
we present our main conclusions. 

We assume a cosmology with Q\ = 0.7, Q,m = 0.3, 
and Ho = 70 km s _1 Mpc - . Magnitudes are in the AB 
system. 

2. OVERVIEW OF METHODOLOGY 

We utilize the COSMOS broad, intermediate, and 
narrow band photometric catalog (jCapak et al.l 120071 : 



I Sanders et al.l 120071 ) for the selection of quasars and 
identification of contaminants. With 29 bands of 
well-matched photometry, these data constitute low- 
resolution spectra for all objects over the wavelength 
range 0.1-8/im (Figure [1]). This is generally sufficient to 
distinguish stars from high-redshift quasars. With spec- 
tra of sufficient quality for all sources in the field, the 
population of quasars can be well-constrained down to 
th e limiting magnitu de. Our approach is similar to that 
of iWolf et ail l|2003j ). who used the 17 filters (5 broad 
and 12 intermediate band) of the COMBO-17 survey to 
identify quasars. COSMOS is of similar resolution to 
the COMBO-17 survey in the optical, but significantly 
deeper and covering a wider wavelength range. 

Because neither the properties of the faint, high- 
redshift quasar population nor the properties of the con- 
taminating stellar population are very well constrained, 
we individually inspect the spectral energy distributions 
and imaging data to identify quasars and reject contam- 
inants, as described in § 5. To ensure that our classifi- 
cations are unbiased, we also cross-check against auto- 
mated x 2 fitting of model SEDs to the low resolution 
spectro-photometry of COSMOS (§ 6.2). The two meth- 
ods produce similar results, but we argue the visual clas- 
sification is more accurate because it allows for a broader 
range of quasar properties as well as the identification of 
photometry errors that can influence the x 2 result. 

The primary method we use can be summarized as 
follows: 

1. Find a coarse selection (based on point-source mor- 
phology, the presence of a Lyman-break, and a 
power law infrared slope) that is highly complete 
for type-1 quasars at z ~ 3 — 5, sacrificing reliability 
to the extent necessary to achieve high complete- 
ness. 

2. Assess the 29-band photometric and imaging data 
for each candidate to remove contaminants. 

3. Cross-check the resulting quasar list with known 
high-redshift quasars in COSMOS to confirm our 
ability to recover quasars with high completeness. 

4. Compute the luminosity function with the resulting 
sample. 

Several factors specific to the COSMOS field make this 
approach reasonable. We list these below. 

• Over 40 confirmed z > 3 quasars are known in 
COSMOS, which were selected in different ways 
(X-ray, infrared, optical). These quasars serve as a 
guide in developing our selection criteria, and pro- 
vide an important check on our ability to distin- 
guish quasars from contaminants using photomet- 
ric and imaging data. 

• A large spectroscopic sample of faint AGN and 
galaxy candidates obtained by Keck and VLT can 
be used to test contamination. 

• High-resolution imaging of the COSMOS field with 
the HST lets us restrict our sample to true point 
sources, limiting the contamination from high- 
redshift star-forming galaxies. 




■ Accurate photomet ric redshifts (jllbert et al.l 120091 : 
ISalvato et aTl r2009) can be used for sources lacking 
spectroscopic confirmation. 

■ Deep Chanda X-ray imaging of the central 0.9 deg 2 
of the COSMOS field allows another check on the 
completeness and reliability of our approach, and 
also permits a comparison of our estimated QLF 
with the X-ray luminosity function at similar red- 
shifts. We present this analysis in § 9. 

3. INITIAL CANDIDATE SELECTION 



^ z — 3.2 quasar SED 

— SDSS quasar template fit 



Wavelength (A) 

Fig. 1. — This confirmed quasar at redshift z = 3.2 illustrates the quality of the low-resolution spectrum afforded by the COSMOS 
photometric data. Filter bandwidth is indicated with horizontal bars. The observed -frame positions of the Ly man break and Lyman limit 
are indicated with vertical blue lines. Overlaid is a typical QSO spectral template (Vandcn Berk et al. 2001). Features of a high-redshift 
quasar SED are clear in the photometry, including the Lyman break at rest-frame 1216A, Lyo emission, broad quasar emission lines, and 
the rising SED into the observed-frame infrared. 

scales with its magnitude. The point source selection is 
shown to be robust to J = 25, which sets the magnitude 
limit for our search. 22383 candidates meet this selection 
criterion. 

Next we use the \ 2 parameter, computed through tem- 
plate fitting in the COSMOS photometric redshift cata- 
log, to separate stars fro m quasars. The Le Phare code 
(jArnouts fc Ilbertl 1201 If) was used to determine photo- 
metric redshifts (see lllbert et al.1 [20091 for details). As a 
result of extensive fitting against templates, sources are 
assigned three \ 2 values: x 2 ga n x 2 so , and xltan whicn are 
computed by fitting galaxy, AGN, and stellar templates, 
respectively, to the source photometry. Because we re- 
strict our search to point sources, a comparison between 
Xg SO and xitar can distinguish most stars from quasars. 
As seen in Figure [H the relaxed cut Xqso/xltar — 2-0 
selects all known z > 3 quasars, while rejecting the ma- 
jority of stars. This leaves us with 8125 candidates. 

The effectiveness of \ 2 fitting is due to two main fac- 
tors. One is the COSMOS intermediate band data, which 
can differentiate sharp features such as the Lyman break 
and emission features in quasar photometry from the 
smooth photometry of stars. The other (more impor- 
tant) factor is Spitzer IRAC photometry. Stars display 
Rayleigh- Jeans black body curves that decline in the in- 
frared, whereas high redshift quasars display power law 
SEDs that are flat or rising in the infrared. 

Next we apply a cut designed to select high-redshift 
quasars. The Lyman break, caused by the scattering of 
radiation at wavelengths shorter than 1216 A by neutral 
hydrogen in the IGM and quasar host galaxy, results in 
a reduction in U band flux for quasars at z > 3.1. In 
the mid-infrared, however, high redshift quasars remain 
relatively bright due to their intrinsic power law spectral 
slopes. Therefore, we accept sources meeting either of 
the following cuts: U — chl > 1.5 or U — ch2 > 1.5. Here 



Spectroscop ic ca mpaigns, e.g. iLillv et al.l (|2007[) . 
iTrump et alT(|2009[ ) , llkeda et al.l (|2011| ). and Capak et al. 
(in preparation) , have confirmed over 40 quasars at z > 3 
in the COSMOS field. These known quasars were origi- 
nally selected in different ways (X-ray, infrared, optical) 
and therefore display a range of SEDs. These confirmed 
sources are used to guide the development of a coarse ini- 
tial selection, which is then tested for completeness with 
simulated quasar photometry, as described in § 4. 

It is worth emphasizing that the selection presented 
here is intended only to be a weak filter against contam- 
ination, in order to reduce the number of candidates to 
a reasonable number while maintaining high complete- 
ness. The careful rejection of contaminants occurs later 
through inspection of the photometric and imaging data, 
as described in § 5. 

We begin by selecting point sources in the magnitude 
range 16 < / < 25, restricting our study to the region of 
COSMOS covered by ACS imaging (1.64 deg 2 ). Spatially 
extended sources as defined in the Hubble Advanced 
Came ra for Surveys (ACS) catalog of iLeauthaud et ahl 
(2007) are excluded. The test for point morphology 
is based on the ratio of the peak surface brightness 
(MILMAX) to the magnitude (MAG_AUTO), and ex- 
ploits the fact that the light distribution of a point source 
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Fig. 2. — Comparison of the minimum \ 2 values for QSO and 
star template fitting to the COSMOS photometry of point sources 
with spectra. The cut shown selects objects with Xqso/Xstar — 2-0, 
which retains all known z > 3 quasars while rejecting a significant 
fraction of stars. The spectroscopically confirmed stars shown here 
are primarily F, M, L, and T dwarfs, which are typical co ntami- 
nants in high-redshift quasar searches ([Richards ct al. 20QU). 

chl and ch2 refer to Spitzer IRAC channel 1 (3.6/xm) and 
channel 2 (4.5/im). These cuts were motivated by the 
spectroscopically confirmed sample as well as simulated 
quasar spectra, and are significantly relaxed in order to 
achieve high completeness. While the cuts are not signif- 
icantly different, we accept sources meeting either one in 
order to avoid losing quasars with bad IRAC photome- 
try in one of the bands. This leaves 6667 candidates. To 
safeguard against losing unusual quasars, we also accept 
sources with photometric redshifts z qso > 3. This adds 
28 candidates, for a total of 6695 candidates after this 
step. 

Finally, we require that the source is detected at a sig- 
nificant level (tuab < 24.0) in either IRAC channel 1 or 
channel 2. Because quasar SEDs typically rise into the 
infrared, this condition will be met by the majority of 
quasars with I < 25, while it serves to eliminate faint 
stars or other unexpected contaminants. The final num- 
ber of candidates after applying this selection is 4009. 

This selection, by utilizing only the Lyman contin- 
uum break and power law infrared emission of high- 
redshift quasars, largely avoids the uncertainties associ- 
ated with broad-band color selection. Color selections in 
the rest-frame optical must correctly account for intrin- 
sic quasar SEDs that can vary substantially in spectral 
slope and broad line emission equivalent widths, the dis- 
tributions of which are not well constrained at high red- 
shift. In addition, this selection should be sensitive to 
broad absorption line (BAL) quasars, which show simi- 
lar mid-infrared to opt ical luminosity ratios as non-BALs 
(|Gallagher et al.ll2007h . 

To summarize, the initial (coarse) selection we use is: 

1. 16 < 7(mag_auto) < 25 

2. Point source based on ACS /-band imaging 

3- Xlso/Xlar < 2-0 



4. U - chl > 1.5 OR U - ch2 > 1.5 OR z qso > 3.0 

5. m AB (3-6nm) < 24.0 OR myuj(4.5/zm) < 24.0 

The 4009 candidates selected include all 39 known z > 
3.1 quasars in COSMOS. 

A comparison of our se lection with the hig h-redshift 
X-ray sources reported in iCivano et al.l (|2011l ) confirms 
that we select all high redshift X-ray sources with / < 25 
and point morphology (27 total). The other 74 sources 
in the X-ray z > 3 sample are primarily classified as ex- 
tended (90%), or lie outside our /-band magnitude limit, 
or both. Keck DEIMOS spectra (Capak et al., in prepa- 
ration) of 12 of the extended z > 3 X-ray sources con- 
firms that these are obscured, type-2 AGN without broad 
emission lines, and therefore not appropriate to include 
in this work. 

These pieces of evidence indicate that the selection out- 
lined above is highly complete for unobscured quasars at 
z > 3.1. This is verified using simulated quasar photom- 
etry, as described in the next section. 

4. SELECTION COMPLETENESS 

We follow the standard procedure of testing the se- 
lection by simulating a large number of quasars in red- 
shift/magnitude space. We use the three composite spec- 
tral templates QSOl, TQ SOl, and BQSOl from the 
SWIRE template library (jPolletta et al.ll2007l) . each rep- 
resenting a different infrared/optical flux ratio, and ex- 
tend them into the far-UV using the composite HST 
spectral template of iTelfer et al.l (|2002h . By letting the 
intrinsic spectral slopes for A > 1216 A and A < 1216 A 
vary randomly acc o rding to th e distributions giv en in 
iVanden Berk et al.l (pOOl and ITelfer et al.l (pOOl . we 
generate 50 QSOl templates and 25 each of the TQSOl 
and BQSOl templates. These 100 base templates are 
checked against the photometry of known quasars in 
SDSS, including IRAC p hotometry for quas ars in the 
Spitzer Heritage Archive (Tep litz et al.l [2012). to verify 
that they are representative of the underlying variation 
in the population of type-1 quasars. 

We create a grid in redshift and / magnitude, with the 
redshift extending from 3 < z < 5 in steps of Az = 0.1 
and the / band magnitude extending from 20 < / < 26 
in steps of Al = 0.1. At each point in this grid we 
generate 20 realizations of each of the 100 model quasar 
templates, adding Poisson and background noise to the 
measured flux in each filter, and c orrecting for in tergalac- 
tic extinction using the model of Madau (1995]). The se- 
lection criteria (4) and (5) are then applied to each of 
the 2000 resulting SEDs to determine the completeness 
at that grid point. 

The result is shown in Figure[3] For / < 23.8, the com- 
pleteness is >90% over the redshift range 3.1 < z < 3.3 
and ^100% at higher redshifts. The marginal incom- 
pleteness for 3.1 < z < 3.3 will be improved by the addi- 
tional inclusion of sources with z qso > 3 in criterion (4). 
There is some incompleteness for 23.8 < / < 25 due to 
criterion (5), which is accounted for when computing the 
luminosity function. 

Figure [3] demonstrates that the photometric criteria 
(4) and (5) are highly complete, but it does not account 
for possible incompleteness introduced by criteria (2) and 
(3). We briefly discuss the effect these may have. 
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Fig. 3. — Completeness map for our photometric selection. For 
magnitudes brighter than / ~ 24, the selection is over 90% com- 
plete over the redshift range 3.1 < 2 < 3.3 and ^100% complete at 
higher redshifts. Because we also accept all sources with z qso > 3, 
the completeness from 3.1 < z < 3.3 will be enhanced over what 
is shown. At /-band magnitudes fainter than / ~ 24 there is some 
incompleteness due to criterion (5), which is accounted for in com- 
puting the luminosity function. 

The third selection criterion (x^so/Xstar ^ 2.0) is based 
on the results of extensive fitting of the COSMOS sources 
against spectral templates. Checking this criterion with 
our modeled quasar templates would not be very reveal- 
ing, as it would constitute checking simulated photome- 
try against simulated photometry. However, the evidence 
in Figure [2] gives us confidence that this very relaxed cri- 
terion introduces negligible incompleteness. 

The point source criterion can potentially induce in- 
completeness because ACS imaging, with a point-spread 
function (PSF) with FWHM of 0.12" at the detector, re- 
solves all but the most compact of objects. While this is 
desirable in that it excludes galaxies from consideration, 
it may also reject intrinsically faint quasars with host 
galaxy light contributing to the overall flux. We discuss 
this issue in more detail in § 8.4. 

5. EXAMINATION OF CANDIDATES AND 
IDENTIFICATION OF QUASARS 

The selection described above is highly complete but 
also unreliable, in the sense that a small fraction of the 
selected objects actually are quasars at the redshifts of 
interest. However, careful examination of the COSMOS 
photometric and imaging data for each of the 4009 can- 
didates is relatively fast and allows us to confidently re- 
ject typical contaminants such as dwarf stars and lower- 
redshift quasars. 

Visual examination of the data allows more robust clas- 
sification than a reliance on SED template fitting. The 
29 bands of photometry and corresponding imaging data 
reveal stars with unusual colors that mimic quasar SEDs 
closely, often because of photometric contamination, as 
well as true quasars that automated tests fail to find be- 
cause of contaminated or unusual photometry In addi- 
tion, coarse redshift estimates can be made based on the 
position of Lya emission and/or the Lyman break, which 
can then be compared with the results from automated 
fitting. 

The relevant high redshift quasar characteristics we 
look for in the COSMOS photometric and imaging data 



include: (1) a strong spectral break at rest-frame 1216 A, 
(2) a power law SED that is flat or rising in the infrared 
IRAC bands, (3) Lya emission evident in intermediate 
band filters, and (4) often other strong rest-frame UV 
emission lines (e.g. CIV A1549, CIII] A1909) evident in 
intermediate band filters as well. 

Stars, in contrast, have SEDs with characteristic 
Rayleigh- Jeans shapes that decline in the IRAC bands. 
The decline in the blue side of a stellar SED is less abrupt 
than the Lyman break in the SED of a high redshift 
quasar, and there are no emission lines in the interme- 
diate bands. Figure |4] illustrates, for some of the most 
difficult cases, the clear differences in the SEDs of high 
redshift quasars and stellar contaminants. 

The candidates wer e examined using Specpro 
([Masters fc Capakl 1201 lh . an interactive IDL tool de- 
signed primarily for spectral analysis in the context of 
multiwavelength surveytQ. This program displays the 
imaging and photometric data for each candidate and 
allows us to perform rough SED fitting and redshift 
estimation. The first author went through all candidates 
after significant calibration against the SEDs of known 
stars and high-redshift quasars, and subsets of the data 
were examined by coauthors to establish the consistency 
and reproducibility of classification. It was found that 
classifications were highly consistent among different 
classifiers. 

As expected, our sample includes a large fraction of 
stellar contaminants. Most of these are dwarf stars, while 
a fraction are stars whose optical fluxes are incorrectly 
matched to a nearby IR-bright source, which is appar- 
ent through examination of the stamp images. The mis- 
match artificially boosts the SED in the IRAC bands, 
such that the Xqso vs - X^tar comparison does not favor a 
star assignment. 

We also find a number of quasars or compact galaxies 
at z < 3.1 in our sample, which reflects the fact that 
our Lyman-break criterion (4) is very relaxed. This can 
also be seen in the slow decline in completeness toward 
z — 3 in Figure |3l Lower-redshift objects are identified 
by the position (or absence) of the Lyman break / Lya 
emission in the SED. Of the contaminants, we estimate 
that 28% are z < 3.1 objects (quasars or compact low- 
z galaxies), and 64% are stars. The remaining 8% are 
unclear based on the photometry but likely fall in one 
of these categories as they do not show characteristics of 
z > 3.1 quasar SEDs. 

For the sources identified as likely high redshift 
quasars, we assign a confidence flag ranging from 2 to 4, 
where 4 is highly confident and 2 is somewhat confident. 
We emphasize that sources even marginally consistent 
with being high-redshift quasars are retained (usually at 
confidence 2) in order to maintain high completeness. 
Candidates assigned a confidence of 4 typically are X-ray 
detected and / or show strong emission lines in intermedi- 
ate bands, as well as a clearly defined Lyman break and 
rising infrared SED. Candidates assigned a confidence 3 
show a well-defined Lyman break and/or Lya emission 
evident in an intermediate band filter as well as a rising 
infrared SED, but are not clearly detected in X-ray. We 
are confident that these are high redshift quasars based 
on their point-source morphology and SED. Candidates 

1 specpro.caltech.edu 
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Fig. 4. — SEDS from the quasar candidates selected in 111. Stars are shown on the left and confirmed quasars on the right. The stars are 
dwarfs represe ntative of those typically contaminating quasar searches, and the quasars are among the faintest known. Stellar templates 
llPicklesI 119981 ) are overlaid in green. The obvious differences between the SEDs of the objects make it relatively simple to distinguish 
them, even without high-resolution spectroscopy. The observed-frame position of the the Lyman break and Lyman limit are indicated on 
the quasar SEDs with vertical blue lines. The distinct break, with a power law SED in the infrared bands, is characteristic of faint, high 
redshift quasars and distinguishes them from stars. 



assigned a confidence 2 show a moderately convincing 
photometric break and are flat or rising in the infrared, 
but do not necessarily show emission features in interme- 
diate bands, nor are they X-ray detected. Of our candi- 
dates these are the most likely to be either very compact 
galaxies or unusual stars rather than quasars. Follow-up 
spectroscopy on these sources may be required to refine 
the results presented here. 

6. QUASAR SAMPLE 

Of 4009 candidates, we find 155 likely type-1 quasars at 
z > 3.1. Of these, 48 are considered confidence 2 sources, 
54 are considered confidence 3, and 53 are considered 
confidence 4. The final list of the 155 confirmed and 
likely z > 3.1 quasars is given in Table|3]in the Appendix. 



6.1. Identification Reliability 

An important point is that the quasar candidates were 
examined without knowing which were the previously 
confirmed quasars, yet every one of the 39 known quasars 
at z > 3.1 in the COSMOS field was correctly classified 
based on the imaging and photometric data. The method 
employed thus independently recovered every previously 
known type-1 quasar in the field. Figure [5] shows that the 
confirmed sample includes a nu mber of very faint s ources, 
including those discovered by llkeda et al.1 (|2011h . This 
is strong evidence that we do not miss actual high- 
rcdshift quasars in the initial sample. If anything, we 
may overcount by misidentifying either stars or high- 
rcdshift galaxies as quasars. Because stars display easily 
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recognizable SEDs, we consider it unlikely that they con- 
stitute a significant contaminant. 

While stars are relatively easy to reject, high redshift 
galaxies display a similar overall SED to a quasar and 
could contaminate the sample. We expect that few of 
the candidates identified as quasars are actually com- 
pact galaxies, because the strict point-source criterion 
will eliminate nearly all galaxies at these redshifts. We 
discuss this potential source of error in more detail in 
§ 8.3. 

6.2. Comparison with Automated x 2 Selection 

As another check on the reliability of our classifica- 
tions, we compare our final sample with a purely au- 
tomated selection using the results of the extensive x 2 
fitting against ste llar and AGN templates described in 
lllbert et all (|2009t ). Of the initial sample of 4009 candi- 
dates, 195 objects have photometric redshift estimates 
above z = 3.1. Limiting the sample to objects with 
Xqso/x 2 tar ^ 1-0 (i- e -j quasar is the favored solution from 
X 2 fitting) the final automated sample is 143 objects, in 
very close agreement with our final sample of 155 from 
visual classification. 

Figure |6] shows that, of the 195 objects with photomet- 
ric redshifts above 3.1, our visual classification typically 
retains objects with Xqso/xltar ^ 1 an d rejects objects 
with x\solx 2 tar ^ 1> as would be expected if the visual 
classification is in agreement with x 2 fitting. 

Of the 143 objects that would be selected automati- 
cally, 91 are in common with our final selection. The 
sources that we rejected but automated fitting would re- 
tain are, for the most part, stars mismatched to a nearby 
IR-bright source, artificially elevating their IRAC fluxes 
and thus lowering their Xqso value. This is reflected in 
the fact that these sources have higher average IRAC-to- 
optical match distances. 

On the other hand, we identify 64 likely quasars that 
are rejected by the automated selection. These primarily 
are rejected due to a low photometric redshift estimate. 
However, they are consistent with being high-z quasars 
based on a weak photometric break and power law in- 
frared spectrum, so we retain them. 

Despite the slightly different samples, it is clear that 
automated selection based on x 2 fitting produces essen- 
tially the same answer as careful rejection of contami- 
nants on an object-by-object basis. 

6.3. X-ray Properties 

The Chandra COSMOS survey (C-COSMOS, 
lElvis et alJ l2009) is a 1.8 Ms Chandra program covering 
the central 0.5 deg 2 of the COSMOS field with an 
effective exposure of ^160 ks, and an outer 0.4 deg 2 
area with an effective exposure of ^80 ks. In addition, 
the entire COSMOS field has been observe d in the X-ray 
with XMM- Newton (|Hasinger et all 120071 ) to a brighter 
X-ray flux limit. 

Of the final sample of confirmed and likely quasars 
presented here, 27 are X-ray detected by Chandra 
(jCivano et al.ll201lL 2012 sub mitted) and 12 are detected 
bv XMM (jBrusa et al.ll2010D . The remainder of the ob- 
jects are undetected and therefore have X-ray fluxes be- 
low the detection limits of either Chandra deep, Chan- 
dra shallow, or XMM, depending on their location in 
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Fig. 5. — The distribution in 7(mag_auto) for the spectroscop- 
ically confirmed subset of our final quasar candidates. We have 
recovered both relatively bright and very faint known quasars in 
the field, indicating that there is high overall completeness in our 
classification method regardless of quasar luminosity. 




Fig. 6. — The fraction of the z qS o > 3.1 candidates we identify 
as high redshift quasars as a function of the ratio x\solx\tar- As 
expected, the lower this ratio, the more likely that we identify 
a candidate as a quasar through inspection of the imaging and 
photometry, showing that our visual inspection is in accord with 
the results of x 2 fitting against spectral templates. 

the COSMOS field. We perform a stacking analysis of 
47 unconfirmed sources without X-ray detection in the 
Chandra region to investigate whether their X-ray fluxes 
are consistent with e xpected quasar em ission. 

We use CSTACK (jMivaii et al.|[2008l ) to compute the 
number of X-ray counts at the optical position of the 
undetected sources. This program correctly accounts for 
the different exposure times depending on source posi- 
tion. We find a detection in the soft-band 0.5-2.0 keV 
(Figure [7]), corresponding closely to rest-frame 2-10 keV 
for our sources, with a significance of ~5-6 sigma, as de- 
termined through random stacking in blank fields. We 
find no clear detection in the observed-frame 2-8 keV 
hard band. The stacked detection in the soft band pro- 
vides evidence that these sources are quasars with indi- 
vidual X-ray fluxes below the Chandra detection limit, 
while the lack of a clear hard band detection suggests no 
significant X-ray obscuration. 

As another check on whether the X-ray flux limits on 
the nondetected candidates are consistent with these ob- 
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Fig. 7. — Stacking analysis of 47 likely quasars that fall in the 
Chandra region but are not individually detected in X-ray. The im- 
ages are 30" X 30" and are smoothed with a Gaussian kernel with 
FWHM of 6". The ~5-6 sigma detection in the soft band (corre- 
sponding to the hard X-ray band in the rest frame at z ~ 3 — 4) 
is evidence that at least a fraction of the individually undetected 
sources arc quasars. The lack of a clear detection in the observed- 
frame 2-8 keV stack indicates they are likely unobscured, intrinsi- 
cally faint objects that fall below the detection limit in the hard 
band. 

jects being quasars, we compute the a ox ratio for sources 
in the C-COSMOS region, making use of the SEDs to de- 
rive the rest-frame UV fluxes. The a ox parameter repre- 
sents the spectral slope between the UV and X-ray flux : 

_ log(/x//uv) m 

Vox = — t -, r (1) 

log(l/x/Z/ uv ) 

where / uv is the 2500 A flux and fx is the 2 keV X-ray 
flux. For X-ray detected sources we derived the 2k eV X- 
ray flux from the 0.5-2 keV flux (|Elvis et al.l 12009ft . For 
sources without X-ray detection, we compute an upper 
limit on the X-ray flux, taking into account the X-ray 
flux limit at the optical position of each source. The re- 
sult in Figure E] shows that the nondetected sources tend 
to have lower UV luminosity, with upper limits on the 
a ox ratio that are generally consistent with previously 
measured correlations between a ox and the 25 00 A lu- 
minosity USteffen et alj|2006l fYoung et al.ll2010f ). In con- 
junction with the stacking results, this gives additional 
evidence that the sources are faint, type-1 quasars with 
X-ray fluxes below the Chandra detection limit. 

7. LUMINOSITY FUNCTION 

We derive the luminosity function for the redshift bins 
3.1 < z < 3.5 an d 3.5 < z < 5 .0, using the standard 
l/V max estimator (|SchmidtHl96l ). 

The luminosity function for quasars is generally found 
to be well-described by a double power law of the form 



<£(M* ) 

lQ0.4(a + l)(M-M* )_|_l()0.4(/3 + l) (M-M* ) 



(2) 



where $(M*) is the normalization, M* 450 is the break 
luminosity between the bright and faint en ds, a is the 
brigh t-end slope and j3 is the faint-end slope (| Boyle et al.l 
1988). We use this parameterization to fit our results 
in con junction with bright-end results of Richards et al. 
(2006) derived in the SDSS. Before presenting the results 
we first outline the determination of redshift and absolute 
magnitude for the sample. 

7.1. Redshift Determination 

Only 39 of the 155 sources in the final sample have 
spectroscopic redshifts. For those that do not, we use 
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Fig. 8. — The X-ray / UV correlation for the sources in our sam- 
ple in the Chandra region of COSMOS. The a ox ratios (or upper 
limits on the ratio for sources without individual X-ray detection) 
are generally consistent with measured correlations at lower red- 
shifts. The dotted lines indicate the 1-sigma dispersion in the a x 
relation of Young et al. 2010. 

the photometric redshift z qso determined by fitti ng AGN 
templates to the source SED (lllbert et al.ll200ll . In the 
cases in which a candidate is X-ray detected but not 
spectroscop ically confirm e d, we adop t the photometric 
redshifts of iSalvato et al.l (|2009l 1201 If ), which have been 
computed for X-ray sources in COSMOS and are opti- 
mized for AGN. 

We check to verify that the photometric redshift is in 
rough agreement with the redshift estimated from visual 
examination of the source SED. If there is a very large 
disagreement between the photometric and visual red- 
shift estimates, we adopt the visually estimated redshift. 
This accounts for occasional failures of the photometric 
redshift determination due to contaminated photometry 
in a small number of bands. 

In Figure |9] we show both the photometric redshift and 
estimated redshift compared with the spectroscopic red- 
shift for the subset of spectroscopically confirmed candi- 
dates. The photometric redshift is generally more accu- 
rate, but contains a small number of outliers. The vi- 
sually estimated redshifts are low on average, but there 
are no outliers. These estimated redshifts were made 
by identifying the position of the Lyman-break/Lyman- 
limit in the source SED. The fact that they are system- 
atically low by a small amount reflects a small bias in 
this visual estimate. We only adopt these values in cases 
in which the photometric redshift is likely incorrect, pri- 
marily as a result of contaminated photometry (29 of the 
80 sources in the 3.1 < z < 3.5 bin and 9 of the 48 sources 
in the 3.5 < z < 5 bin). Because our redshift bins are 
relatively large, the errors introduced by doing so will 
not have a significant effect on the derived luminosity 
functions. 

7.2. Absolute Magnitude Determination 

We follow the convention of giving the quasar lumi- 
nosity function in terms of the absolute magnitude at 
rest- frame 1450 A. To estimate M1450, we fit quasar tern- 
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Fig. 9. — Top: the photometric versus spectroscopic redshift for 
the sources in our sample with spectroscopic confirmation. The 
scatter is ~0.7, with a few outliers. Bottom: the estimated red- 
shifts versus spectroscopic redshifts for the same sources. Our vi- 
sual estimates are systematically low, but with no outliers. These 
values are only adopted when no spectroscopic redshift is available 
and the photometric redshift is clearly wrong. This amounts to 29 
out of 80 sources in the 3.1 < z < 3.5 bin and 9 out of 46 sources 
in the 3.5 < 2 < 5 bin. 

plates to the source SED, which accounts for elevated 
fluxes due to quasar emission lines. We found that fit- 
ting to the measured magnitudes in the three broad-band 
Subaru filters V, R, and I and interpolating to find the 
observed-frame 1450 A magnitude often overestimates 
the luminosity due to the contribution of quasar emis- 
sion lines to the measured broad-band magnitudes. 

When performing the SED fits to find M1450, only fil- 
ters within ±2000 A of the observed-frame wavelength 
corresponding to rest-frame 1450 A are considered. The 
-M1450 values we find through SED fitting are typically 
~0.2 magnitudes fainter than what we found by extrap- 
olating the broad band magnitudes, but with significant 
scatter (Figure ITU1) . 

Given the observed magnitude m corresponding to 
rest-frame 1450 A radiation, the absolute magnitude is 
found with 

M 1450 = m - 51og(di/10) + 2.51og(l + z) (3) 

The second term on the right-hand side is the distance 
modulus; the final term is a correction for the effective 



TABLE 1 

Binned luminosity function at 2 ~ 3.2 and 2 ~ 4. 

AMi4so N QSO * (xl0~ 7 ) A$ (xlO" 7 ) 
(Mpc -3 mag -1 ) 

z ~ 3.2 



[-25.5,-23.5] 
[-23.5,-22.5] 
[-22.5,-21.5] 
[-21.5,-20.5] 

Z r~J 4 



4 
9 

26 
11 



5.3 
11.9 
39.0 
59.8 



2.6 
3.9 
7.1 
8.4 




5 
9 

14 
20 



1.9 
3.5 
6.4 
20.6 



0.8 
1.2 
1.4 
2.6 



narrowing of the frequency interval sampled in the ob- 
served frame. 
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Fig. 10. — Comparison of the M1450 value derived from a linear 
fit to broad band filters vs. that derived from SED fitting. The 
absolute magnitude estimated using the SED fit is fainter by ~0.2 
magnitudes on average, as indicated with the dotted line. This 
reflects the fact that the broad band fluxes can be elevated relative 
to the underlying continuum by broad emission lines. 



7.3. The Luminosity Function at z ~ 3.2 and z ~ 4 

We compute the luminosity fuction using the derived 
redshifts and absolute magnitudes for the sample, cor- 
recting for the incompleteness at / > 24 from the selec- 
tion function in Figure [3l The luminosity function we 
derive for the redshift bin 3.1 < z < 3.5 is shown in 
Figure fTTl There is good continuity with the bright-end 
result derived in the SDSS, as well as close agreement 
with prior result s that have probe d the faint end at simi- 
lar redshift (e.g. iSiana et al.l2008llBongiorno et~aT]|2007l 



IWolf et al.1 12003'). 
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We fit the z ~3.2 QLF to the double power law rela- 
tionship (Equation 2) using our faint-end result in con - 
junction with the bright end from lRichards et al.l ([2006) . 
The best-fit parameters we derive are <£* = 2.65(±2.22) x 
10" 7 , Mf 450 = -25.54 ±0.68, a = -2.98 ± 0.21, and 
13 = -1.73 ±0.11. 

iir 4 ^ — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — ^ 

; + This Work 

1 SDSS (Richards et al. 2006) 
1 SWIRE (Siana et al. 2008) 




Fig. 11. — The z ~ 3.2 luminosity function derived here compared 
with previous results. The solid black line is our best fit to the 
double power law parameterization including the bright-end result 
from the SDSS. Our result is in close agreement with previous 
determinations of the QLF at similar redshift. The faint-end slope 
we derive is ft = —1.73 ± 0.11. 

In Figure [T2] we show the result for the z ~ 4 luminos- 
ity function, together with the previous results from Gil 
and 111. The best-fit parameters to the double power law 
paramaterization are $* = 7.5(±23.5) x 10~ 8 , M* 450 = 
-25.64 ± 2.99, a = -2.60 ± 0.63, and /3 = -1.72 ± 0.28. 
Our faint end result is in agreement with the result re- 
ported by 111, and the faint-end slopes we find at both 
z ~ 3.2 and z ~ 4 agree with th e faint-end s l ope a t 
2 < z < 3.6 (f3 = -1.70) derived in I Jiang et ail (|2006|) . 
However, the space density of faint quasars at z ~ 4 we 
find is lower than the result reported in Gil by a factor of 
~3-4. The source of this discrepancy is unclear, but may 
be due to contamination (dwarf stars and high-redshift 
galaxies) in the sample of Gil at the faintest magni- 
tudes probed in that work. The point-source selection 
used in Gil is done with ground-based imaging, which 
is not as effective as HST imaging at resolving compact, 
high-z galaxies. Additionally, the number of stellar con- 
taminants in broad band quasar searches grows rapidly 
at the faintest magnitudes, where spectroscopic follow- 
up is the most challenging and thus the contamination is 
not as well-constrained. 

In Figure [T31 we illustrate the strong evolution of the 
luminosity function from z ~ 3.2 to z ~ 4. We find that 
the space density of faint quasars decreases by a factor of 
four between z ~ 3.2 and z ~ 4, in close analogy with the 
evolution of brighter optical and X-ray quasars at these 




Fig. 12. — The z ~ 4 luminosity function derived here compared 
with previous results. The solid black line is our best fit to the 
double power law parameterization including the bright-end result 
from the SDSS. The dotted green line is the best fit from 111 and 
the dashed red line is the best fit from Gil. The disagreement 
between the bright-end result of Gil and the SDSS data points 
arises because Gil recompute the SDSS results after finding the 
absolute magnitudes directly from the SDSS spectra. Our result is 
in close agreement with that of III. The disagreement in the faint- 
end space density between the result derived here and that derived 
in the DLS+NDWFS is significant, with that result higher by a 
factor of ~4. 



TABLE 2 

Double power law parameters, combining SDSS and this 

WORK. 



*(M*) Af* 450 a p 



z ~ 3.2: 

2.65(±2.22) x 10~ 7 -25.54 ±0.68 -2.98 ± 0.21 -1.73 ±0.11 
z ~ 4 : 

7.5(±23.5) x 10" 8 -25.64 ±2.99 -2.60 ± 0.63 -1.72 ±0.28 



redshifts. Th is trend is a l so veri fied by the recent results 
presented in llkeda et all (|2012f ) , where it is found that 
the space density of faint quasars at z ~ 5 is significantly 
lower than at z ~ 4. 

8. POSSIBLE SOURCES OF ERROR 
8.1. Cosmic Variance 

The comoving volume probed here at z ~ 4 is suf- 
ficiently large (2.6 x 10 7 Mpc 3 ) that we expect cosmic 
variance to ha ve little effect. Ho wever, it has been shown 
in the SDSS (|Shen et all 120071 ) that the bias of lumi- 
nous, high redshift quasars relative to dark matter halos 
increases dramatically with redshift, so we attempt to 
estimate the likely error introduced by cosmic variance. 

We make use of the cosmic variance cookbook of 
iMoster et al.l (|2011l ). in which the bias of galaxies as a 
function of stellar mass and redshift is used to compute 
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Fig. 13. — The evolution of the luminosity function from z~4 
to z~3.2, combining the bright end derived from the SDSS 
(Richards ct al. 2006) with the faint end derived here. The space 
density of faint quasars increases by a factor of ~4 from z~4 to 
2~3.2. Overplotted are the best-fit lines for each luminosity func- 
tion using the double power law parameterization. The faint-end 
slopes we derive for z~3.2 and z~4 are /3 = —1.73 ± 0.11 and 
P = -1.72 ± 0.28, respectively. 

the counting error for galaxies due to cosmic variance. 
Quasars are expected to be significantly biased with re- 
spect to the matter distribution, existing preferentially 
in the highest mass halos. This leads us to make the ap- 
proximation that quasars are similarly biased to the most 
massive (M > 10 11 M Q ) galaxies. The bias factor esti- 
mated in this way is i n relatively good agreement with 
the results presented in lShen et al.l (|2007l) for quasar clus- 
tering in the SDSS (and may actually be slightly overes- 
timated). 

In the redshift window of 3.5 < z < 5 in the COS- 
MOS field, we compute an uncertainty of ~19% in our 
quasar count, indicating that cosmic variance will con- 
tribute only modestly to the derived luminosity function. 

8.2. Contamination by Stars 

The consistent, well-matched wavelength coverage of 
the COSMOS data, extending from the UV to the in- 
frared, enables the effective removal of stellar contami- 
nants, the SEDs of which differ substantially from those 
of quasars. Figure [6] illustrates that our final selection 
preferentially retains objects for which \ 2 fitting disfa- 
vors stellar templates, and the close agreement with prior 
studies at z ~ 3 also indicates that contamination is low. 
We note that, even it there were significant stellar con- 
tamination in the final sample, this would imply that the 
luminosity functions are upper bounds, making it more 
difficult to reconcile the result at z ~ 4 found here with 
that of Gil. 

8.3. Contamination by Star-forming Galaxies 

Because we only have spectroscopic confirmation for a 
fraction of our sample, we may overestimate the faintest 



end of the luminosity function by inadvertantly counting 
compact, star-forming galaxies as quasars. To estimate 
the contamination from high-redshift galaxies, we ana- 
lyze a sample of 386 confirmed star-forming galaxies at 
z > 3 in COSMOS. These sources were selected as part of 
a large spectroscopic campaign with Keck DEIMOS (Ca- 
pak et al., in preparation) aimed at finding high-redshift 
galaxies down to I = 25 from 3 < z < 6. A number 
of techniques (Lyman-break criteria, narrow-band selec- 
tion, etc.) were employed to achieve high completeness, 
with no selection based on morphology. By applying our 
quasar selection to these objects, we can get a reasonable 
sense of the contamination due to compact galaxies. 

Of the 386 confirmed high-redshift galaxies, only 2 
(0.5%) pass our selection (and were removed from our 
final sample), with the rest rejected primarily on the 
basis of morphology. We can combine the fraction of 
high-redshift star forming galaxies likely to contaminate 
our sampl e with the UV lumino sity function for galaxies 
at z ~ 4 ()Bouwens et al.ll2007t ) to estimate the level of 
contamination. Restricting our attention to the faintest 
absolute magnitude bin (M1450 = —21) examined at 
z ~ 4, we estimate a contamination level of roughly 
10~ 6 (Mpc~ 3 mag^ 1 ) by star-forming galaxies. This is 
approximately half the space density we derive in this 
magnitude bin. Therefore, this faintest data point on 
the luminosity function may be elevated somewhat sig- 
nificantly by star-forming galaxy contamination, leading 
to an overestimation of the faint-end slope. This issue 
may need to be settled by follow-up spectroscopy on the 
faintest sources in the sample. 

8.4. Faint Quasars in Bright Host Galaxies 

Conversely, we may underestimate the faintest bins 
of the luminosity functions by rejecting bright galaxies 
hosting faint quasars. This will be the case if a signifi- 
cant fraction of faint type-1 quasars (M1450 > —23) exist 
in galaxies bright enough in the rest-frame UV to make 
the source appear extended in ACS /-band imaging. Two 
pieces of evidence lead us to believe that this is not the 
case. 

First, we examine 12 Keck DEIMOS spectra of Chan- 
dra sources at high redshift that are classified as extended 
in the ACS catalog, and find that these sources do not 
show broad lines indicative of a type-1 quasar. Therefore 
there is no evidence from spectroscopy of high-redshift 
X-ray sources of a population of type-1 quasars that we 
miss by restricting our search to point sources. On the 
contrary, high-redshift X-ray sources that are morpho- 
logically extended seem to be exclusively type-2 AGN. 

Second, we note that, in the sample of ^1400 high- 
redshift candidate galaxies in COSMOS followed up on 
with Keck DEIMOS spectroscopy, we do not find cases 
of broad-line objects with extended morphology. If a 
significant number of faint type-1 quasars at z ~ 4 have 
host galaxies bright enough in the rest-frame UV to make 
the galaxy appear extended in ACS imaging, we should 
find serendipitous examples of these in a large sample of 
candidate high-redshift galaxies. In fact we do not; high- 
redshift, broad- line emitting objects seem invariably to 
be point sources in the rest-frame UV. 

We conclude from these two tests that we do not miss a 
substantial fraction of type-1 quasars due to our morpho- 
logical criterion, even at the faintest quasar magnitudes. 
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9. COMPARISON OF THE X-RAY AND UV QUASAR 
LUMINOSITY FUNCTIONS AT HIGH REDSHIFT 

The X-ray luminosity function at z > 3 in COSMOS 
was determined re cently for the 0.9 d eg 2 Chandra re- 
gion in COSMOS (|Civano et all (|20Tfl) . hereafter Cll). 
There it is shown that the population of X-ray luminous 
AGN decreases rapid ly above z = 3, in agreement with 
previ ous results (e.g. iBrusa et al.ll2009t iSilverman et al.l 
2008|). In fact, Fi gure 4 of Cll shows a decline in space 
density of X-ray luminous (L x > 10 44 ergs s _1 ) AGN by 
a factor of ~4 between z ~ 3.2 and z ~ 4, in agreement 
with the trend we have demonstrated here for the faint 
end of the UV QLF at the same redshifts. 
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Fig. 14. — The UV luminosity function results, converted to X- 
ray (black diamonds) , compared to the X-ray luminosity fu n ctions 
for X-ray bright Chandra sources derived in [Civano et ah l|2011fl 
including uncertainties (red circles and y ellow shaded area ) and 
for XMM-COSMOS sources (green circles, [Brusa et al.ll2"009f) . The 
space density is roughly a fourth of that found in X-ray, implying 
that ~75% of bright quasars are optically obscured at these red- 
shifts . The blue lines represent the prediction from the G illi et al.l 
(2007) model for obscured and unobscured sources (solid) and un- 
obscured sources only (dashed). 

A rough sense of the obscured fraction at these red- 
shifts can be obtained simply by comparing the high- 
rcdshift X-ray sample with our list of type-1 quasars. 
Doing so, we find that our selection finds 27/101 sources 
in the high-redshift X-ray sample. The X-ray sources 
that we miss are primarily classified as extended in the 
ACS /-band imaging, or are fainter than / = 25, or both. 
As mentioned previously, an examination of 12 spectra of 
high-redshift X-ray sources that are classified as extended 
in ACS imaging shows that they are consistent with be- 
ing type- 2 AGN rather than unobscured quasars. As- 
suming that the X-ray sources our selection misses are in 
fact optically obscured, this leads to an estimated type-2 
fraction of roughly 75%. 

To more quantitatitvely compare our results with those 
found in Cll, we convert our rest- frame UV luminosity 
functions to X-ray using the a ox relationship presented 



in ISteffen et all (120061) : 

log(Z 2 kev) = (0.721 ± 0.011)log(Z 



(4.531 ±0.688) (4) 



We modify this relationship slightly to put it in terms of 
the luminosity at 1450 A by assuming a typical quasar 
spectral slope. With this relationship, the space density 
of quasars with L x > 10 44 15 erg s _1 is found by inte- 
grating the UV luminosity functions over the absolute 
magnitude interval M1450 = [— oo, —23.45]. The uncer- 
tainty in the conversion is based on the uncertainty in 
the parameters of the a ox relation. 

Our results are plotted as black diamonds in Figure El 
together with the X-ray results of Cll. The same trend 
with redshift is observed, but the points derived from the 
UV QLF are lower by a factor of ^4, in agreement with 
the ratio found by directly comparing the two samples. 
While the Cll space density is in agreement with the 
prediction for obscured and unobs c ured sources (solid 
line) derived from the iGilli et all (|2007| ) X-ray back- 
ground synthesis model, the space density derived from 
the UV QLF is more consistent with the prediction for 
unobscured sources only (dashed line) . Based on the off- 
set, we again find that the fraction of obscured (but not 
Compton-thick) quasars is roughly 75% at z ~ 3 — 5. 

It is interesting to compare this result with th e redshift 
evolut ion of the type-2 fraction determined by lHasingerl 
(2008) (hereafter H08). Figure 9 of that work clearly 
indicates an upward trend in obscured fraction with red- 
shift, which seems to saturate around z = 2 at ^60%. 
This result arises after a correction for the fact that ob- 
servational bias makes it more difficult to identify ob- 
scured AGN at higher redshifts. The type-2 fraction we 
find at z ~ 3 — 4, ^75%, is higher than the result in H08 
and may indicate that the trend of increasing obscured 
fraction with redshift continues to higher redshifts. 

In should also be noted that the result in Figure 9 of 
H08 is actually a normalization of a rather steep gradi- 
ent of obscuration versus X-ray luminosity. The obscured 
fraction decreases with X-ray luminosity at a given red- 
shift, such that for sources above L x = 10 44 at z ~ 3 — 4 
the estimated obscured fraction is actually closer to 20- 
30%, significantly lower than what we find for compara- 
bly luminous sources at these redshifts. Follow-up spec- 
troscopy on more of the z > 3 X-ray sources in COSMOS 
would allow a more firm determination of the type-2 frac- 
tion at these redshifts. 

10. THE IONIZING BACKGROUND DUE TO QUASARS 

Given the luminosity function $(Mi45o, z), we can es- 
timate the rate of production of ionizing photons by 
quasars at redshift z. T o do so, we use result s from 
iVanden Berk et all (pOOll ) and lTelfer et~aTl (pOOl ) to ap- 
proximate the spectral energy distribution of a typical 
quasar: 



L{y) oc 



-°- 5 (1216 A < A < 2500 
- 1 - 76 (A < 1216 A) 



(5) 



The number of hydrogen-ionizing photons emitted per 
second by a quasar of absolute magnitude M1450 is given 

by 

f tyd L(JW). 
«Hl(Ml450) = / 7 dv (6) 

Jlryd hV 
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where the cutoff at 4 ryd is due to pre ferential Hell 
absor ption of photons above this energy (|Madau et al.l 
1999). Similarly, the rate of emission of Hell ionizing 
photons is given by 



L{M 1450 ,is) 



dp 



(7) 



4ryd 



Combining equations (5) and (6) with the QLF, one 
can compute the number of HI and Hell ionizing photons 
produced per second per unit comoving volume at a given 
rcdshift: 



Mm{z)= / $(M,z)h m {M)dM 



A/"HeiiO)= / $(M, z)n HeU (M) dM 



(8) 



(9) 



We use the results for the QLF at z ~ 3.2 and z ~ 4 
presented in § 7 to estimate these ionizing backgrounds 
and their evolution with redshift. 

10.1. HI Ionizing Background 

With the QLFs derived in § 7, we find quasar emission 
rates of hydrogen- ionizing photons at z ~ 3.2 and z ~ 4 
ofTVki = 3.3xl0 50 s" 1 Mpc -3 and 1.0 x 10 50 s" 1 Mpc" 3 , 
respectively. These rates are lower than needed to reion- 
ize the IGM at either redshift. At 2=4, for example, 
the estimated required rate for reionization is A/hi = 
2.4 x 10 51 s" 1 Mpc" 3 (jMadau et alj [l999). more than 20 
times higher than what we find. In addition, the rates we 
compute ignore escape fraction and are therefore upper 
limits on the photons available to ionize the IGM. 

The emissivity of UV radiation due to quasars can 
be found with our lumi nosity functions and com pared 
with the result given in lHaardt fe Madaul (|2012f ). We 
integrate to find the emissivity at 1450 A and use 
a conversion factor of f ° / f « = 0.58 to con- 

912A ' 1450A 

vert to the corresponding emissivity at 912 A. At 
z ~ 3.2 and z ~ 4 we find e « = 4.4 x 10 24 

912A 

" Hz -1 , respectively. 
4.3 x 10 24 and 



and 1.4 x 10 24 e rgs s 1 Mpc 



912A 



lHaardt fc Madaul (120121) find 

2.0 x 10 24 ergs s _1 Mpc -3 Hz -1 at the same redshifts. We 
are in relatively close agreement, although our value at 
z ~ 4 is somewhat lower, likely due to the rapid decline 
in the space density of faint quasars we find. 

The emissivity we find at z ~ 4 is significantly lower 
than that required to maintain reionization at z ~ 7 
(jRobertson et al.l I2010T ) . even assuming a very high es- 
cape fraction for ionizing photons. Therefore, assuming 
that the faint-end of the QLF continues to decline at 
higher redshifts, type-1 quasars can not contribute signif- 
icantly to cosmic reionization of hydrogen at z ~ 6 — 10. 

10.2. Hell Ionizing Background 
We find a QSO emission rate of Hell-ionizing photons 
at z ~ 3.2 of AhoII = 3.1 x 10 49 s" 1 Mpc -3 . At z ~ 4 
we find AhcII = 1.0 x 10 49 s" 1 Mpc" 3 . The significant 
increase in the production by quasars of Hell-ionizing 
(>54.4 eV) photons from z ~ 4 — > 3 is qualitatively con- 
sistent wltli_thfiobs^rved_pj ; tc^ of Hell at 
z ~ 3 (Sh ull et alJlMTAlBolton et al.ll2006l) . While AGN 



are considered the only sources capable of producing suf- 
ficient photons with energies above 54.4 eV to be respon- 
sible for the reionization of Hell, significant uncertainty 
exists with r egard to the speci fics of the Hell reioniza- 
tion process (jShull et al.ll2010T ). This is mainly due to 
degeneracies between temperature, density, and ioniza- 
tion fraction in the IGM, as well as the dependence on 
the spatial distribution and spectral hardness of the ion- 
izing sources. Therefore, we do not attempt to quantify 
the ionizing background needed to produce the observed 
Hell reionization at z ~ 3. 

11. SUMMARY & CONCLUSIONS 

We have searched the COSMOS field for type-1 quasars 
at 3.1 < z < 5 and Jab < 25, with care taken to achieve 
high completeness. Our method has been to apply a weak 
quasar selection that is nearly complete but unreliable, 
and then identify and remove contaminants through care- 
ful inspection of the COSMOS photometric and imaging 
data. This method exploits the fact that stellar con- 
taminants are easily identified with the 29 bands of deep 
photometry spanning the UV to the infrared. We re- 
cover all 39 previously known type-1 quasars at z > 3.1 
in COSMOS as well as 116 additional likely quasars, and 
have presented evidence based on simulations and auto- 
mated classifiers that this sample is highly complete for 
type-1 quasars above z = 3.1 in the HST-ACS region 
(1.64 deg 2 ) of COSMOS. 

Here we summarize our main conclusions. 

1. We find 155 likely type-1 quasars at z > 3.1 in the 
COSMOS field, 39 of which have been previously 
confirmed, down to a limiting magnitude Iab = 
25. We present strong evidence that this quasar 
sample is nearly complete over the redshift range 
3.1< z < 5. 

2. We use the quasar sample to compute the QLF at 
z ~ 3.2 and z ~ 4. The faint-end results we ob- 
tain show continuity with the br ight-end results re- 
ported bv lRichards et al.l (|2006f ) for the SDSS, and 
demonstrate a clear decrease in the space density 
of faint quasars, by roughly a factor of four, from 
2 ~ 3.2 to 2 ~ 4. The faint-end we derive at z ~ 4 
is in good agreement with the result reported for 
the COSMOS field in 111, but disagrees with the 
result reported in Gil. While the source of this 
discrepancy is unclear, we suggest that it may be 
caused by contaminating stars/high- redshift galax- 
ies at the faintest magnitudes in the sample of Gil. 

3. We find no evolution in the faint-end slope /? over 
the redshift range investigated. The faint end 
slopes we find are /3 — —1.73 ± 0.11 at z ~ 3.2 
and (3 = —1.72 ± 0.28 at z ~ 4, similar to those 
found at lower redshift. 

4. We compare the luminosity functions derived here 
with the X-ray luminosity function at z > 3 in 
COSMOS. The optical QLF and the X-ray QLF 
evolve similarly between z ~ 3.2 and z ~ 4 in 
the sense that both show a rapid decline in space 
density. However, the different normalizations of 
the LFs imply that roughly 75% of AGN with 
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log(i x ) > 43.75 are obscured, type-2 quasars at 
z ~ 3 - 4. 



5. We compute the ionizing background due to 
quasars, both of Hell and HI. The rapid increase in 
production of Hell ionizing photons from z ~ 4 — > 
3 that we find is qualitatively consistent with the 
observed onset of Hell reionization at z ~ 3. How- 
ever, given the decline of the faint end of the QLF 
to higher redshift, we conclude that faint quasars 
are unlikely to contribute substantially to cosmic 
reionization of hydrogen. 
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TABLE 3 Summary of the likely quasars found here and used to compute the luminosity function. We give z es t, the redshift we 
estimated based on visual examination, z qao , the photometric redshift from COSMOS, z spec , the spectroscopic redshift, if it exists, and 
z uscd) the redshift we adopt in computing the luminosity function. In addition, we list the confidence flag we assigned the source based 
on visual examination, whether or not the source is detected in X-ray, and the absolute magnitude M1450 we derive. 



ID RA DEC I AB z eat z qso z spec z uaed Flag X-ray Mi 450 

(J2000) (J2000) (AUTO) 

381470 149.85396 1.753672 22.55 3.30 0.11 -99 3.30 4 no -22.11 

1625825 150.43993 2.703496 23.05 3.43 0.05 3.46 3.46 4 yes -21.W 

1199385 150.39026 2.445338 24.07 3.18 2.78 -99 3T8 4 no -20.60 

1575750 150.73715 2.722578 22.23 3.20 3.25 3.32 3.32 4 yes -22.78 

910592 150.73557 2.199578 20.22 3.00 3.45 3.50 3.50 4 yes -25.06 

1163086 150.70377 2.370019 22.44 3.10 3.79 3.75 3.75 4 yes -23.00 

1159815 150.63844 2.391350 22.16 3.00 3.62 3.65 3.65 A yes -22.98 

1605275 150.62006 2.671402 21.83 3.10 3.13 3.14 3AA 4 yes -22.82 

887716 149.49590 1.968019 22.61 3.30 3.23 -99 3.23 4 no -22.38 

1381605 150.57439 2.552482 23.60 3.20 3.13 -99 3.13 4 no -21.54 

503666 150.55476 1.904693 23.90 3.00 3.33 -99 3J33 A no -21.22 

299482 150.49628 1.638000 23.99 3.30 3.35 -99 3J35 A no -21.31 

1856470 150.47568 2.798362 21.28 3.50 3.46 4.12 AA2 A no -24.37 

507779 150.48563 1.871927 22.03 4.30 4.41 4.45 4^45 A no -23.78 

739700 150.45497 1.967424 24.03 3.40 3.48 3.48 3 1 48 A yes -21.20 

504209 150.42476 1.900266 24.79 3.40 1.18 -99 3^40 A yes -20.10 

748601 150.38388 2.074549 23.46 3.50 3.45 -99 3^45 A yes -21.76 

1208399 150.25954 2.376141 21.42 3.30 2.68 3.72 3.72 4 yes -24.47 

1657280 150.24078 2.659058 22.33 3.00 3.33 3.36 3.36 A yes -22.74 

1448618 150.21497 2.582674 2413 5T20 5T26 5T30 OD A yci -22.73 

1463661 150.20885 2.481935 20.11 3.30 3.28 3.33 3.33 4 yes -25.4TT 

1224733 150.20898 2.438466 21.15 3.40 3.50 3.71 3.71 4 yes -24.28 

560579 150.21025 1.853839 22.12 3.40 3.45 -99 3^45 A yes -23.10 

1231613 150.21075 2.391473 22.55 2.90 3.08 3.10 3T0 A yes -22.33 

565710 150.21199 1.818717 24.08 3.50 3.48 -99 3^48 A no -21.18 

348143 150.17534 1.649839 24.19 3.30 3.28 -99 3J28 A no -20.65 

1900976 150.15880 2.808512 23.06 3.50 0.83 -99 3J30 A no -21.58 

1465836 150.13036 2.466012 23.04 3.40 3.45 3.86 3 1 86 A no -22.56 

1221992 150.13364 2.457429 24.14 3.10 3.13 -99 3T3 A yes -20.74 

1226535 150.10098 2.419435 22.33 3.60 4.51 4.66 Am A yes -23.46 

330806 150.10738 1.759201 22.55 3.70 3.92 4.14 4.14 A yes -22.85 

780598 150.09685 2.021495 24.33 3.50 3.23 -99 3J23 A yes -20.49 

582598 150.04268 1.872161 23.82 3.30 3.33 3.36 3.36 A yes -21.80 

585189 150.00940 1.852637 23.67 3.00 3.33 -99 3J33 A yes -21.92 

804307 150.00438 2.038898 21.86 2.80 3.45 3.50 ^50 A yes -23.56 

1261211 149.90547 2.354014 23.19 3.20 3.30 3.27 3J27 A no -21.92 

1054048 149.87920 2.225839 22.70 3.00 3.65 3.65 3.65 A yes -22.77 

1249763 149.89417 2.432972 23.29 3.20 3.35 3.38 3 1 38 A yes -21.84 

1271385 149.86966 2.294046 21.67 3.40 3.41 3.35 3J35 A yes -23.42 

1046585 149.85153 2.276400 23.02 3.30 3.29 3.37 3.37 A yes -22.39 

1511846 149.84576 2.481679 23.11 3.20 3.33 3.36 3.36 A yes -21.94 

1513806 149.78207 2.471342 23.70 3.40 3.26 -99 3 1 26 A yes -20.88 

1719143 149.75539 2.738555 22.87 3.40 3.52 3.52 3.52 A yes -22.26 

1272246 149.78381 2.452135 23.70 4.90 4.95 5.07 5.07 4 yes -22.23 

1273346 149.77692 2.444306 22.78 4.10 3.95 4.16 4T6 A yes -22.65 

1284334 149.77104 2.365819 24.64 3.60 3.45 -99 3A5 A yes -20.25 

1060679 149.73622 2.179933 23.45 3.70 4.24 4.20 4^20 A yes -22.23 

1720201 149.74873 2.732016 23.95 3.10 3.15 -99 3T5 A yes -20.54 

422327 149.70151 1.638375 22.41 3.20 3.17 3.20 3.20 4 no -22.49 

1743444 149.66605 2.740230 22.51 3.20 3.15 -99 3.15 4 no -22.54 

1551171 149.52461 2.531040 22.74 3.30 3.27 -99 3J27 A no -21.71 

1330271 149.52908 2.380164 20.80 3.20 3.09 3.10 3T0 A yes -24.04 

699705 150.58488 2.081361 23.12 3.00 3.29 -99 3J29 A yes -21.83 

1371806 150.59184 2.619375 22.83 3.00 3.12 -99 3.12 3 no -22.10 

1628943 150.55078 2.682909 23.91 3.30 0.45 3.56 3.56 3 no -21.31 

361333 149.92613 1.724807 23.46 3.10 3.41 3.30 3J30 3 no -21.57 

1518518 149.69594 2.603082 24.24 3.30 3.29 -99 3J29 3 no -20.64 

388046 149.72974 1.704067 23.62 3.30 0.17 -99 3J30 3 no -21.32 

1420590 150.36694 2.616221 23.65 3.30 3.83 -99 3J30 3 no -21.37 

990122 150.29726 2.148846 20.46 2.90 3.47 3.33 3J33 3 yes -24.62 

1648871 150.31337 2.716225 24.70 3.70 3.63 -99 3J33 3 no -20.52 

530538 150.27693 1.885094 24.08 4.80 4.87 4.72 4.72 3 no -22.20 

599355 149.85320 1.920494 22.66 3.40 0.74 -99 3A0 3 no -22.17 

1870817 150.24617 2.858176 23.82 3.40 3.55 -99 3.55 3 no -21.45 

1454738 150.20822 2.539911 22.46 3.30 3.24 -99 3.24 3 no -21/5? 

618000 149.83923 1.793837 24.63 3.40 3.55 -99 3J>5 3 no -20.77 

1730531 149.84322 2.659095 22.90 3.20 3.51 -99 3JU 3 no -22.15 

Continued on Next Page. . . 
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TABLE 3 - Continued 



ID 


RA 
(J2000) 


DEC 
(J2000) 


Iab 
(AUTO) 


Zest 


Zqso 


Zspec 


z used 


Flag 


X-ray 


Mi4 50 


1970813 


149.64572 


2.797057 


23.24 


3.40 


4.15 


-99 


3.40 


3 


no 


-21.63 


115356 


150.14893 


1.587745 


24.06 


3.80 


3.88 


-99 


3.88 


3 


no 


-21.50 


1717612 


149.76262 


2.749090 


24.81 


3.90 


0.49 


-99 


3.90 


3 


no 


-20.98 


329051 


150.16891 


1.774590 


23.24 


4.30 


4.35 


-99 


4.35 


3 


no 


-22.71 


790476 


150.11644 


1.963943 


23.90 


3.50 


3.84 


3.41 


3.41 


3 


no 


-21.25 


1326857 


149.55511 


2.402630 


24.32 


3.30 


3.94 


-99 


3.94 


3 


no 


-21.53 


803073 


150.05769 


2.046285 


24.88 


3.16 


0.00 


-99 


3.16 


3 


no 


-21.51 


373838 


150.02069 


1.639405 


24.74 


3.34 


3.20 


-99 


3.20 


3 


no 


-20.43 


374612 


149.91289 


1.632472 


24.07 


3.12 


0.06 


-99 


3.12 


3 


no 


-19.52 


1727644 


149.85581 


2.681164 


24.47 


3.16 


0.12 


-99 


3.16 


3 


no 


-20.71 


376803 


149.85549 


1.768740 


21.06 


3.30 


0.09 


-99 


3.30 


3 


no 


-20.20 


1283260 


149.78825 


2.372171 


23.38 


3.40 


0.24 


-99 


3.40 


3 


no 


-21.61 


594007 


149.85205 


1.959338 


24.39 


3.26 


3.35 


-99 


3.35 


3 


no 


-20.59 


1155791 


150.66118 


2.414451 


24.05 


3.22 


3.15 


-99 


3.15 


3 


no 


-20.97 


1153009 


150.58006 


2.439344 


24.23 


3.16 


3.10 


-99 


3.10 


3 


no 


-20.78 


842269 


149.58212 


2.117000 


23.36 


3.16 


0.11 


-99 


3.16 


3 


no 


-21.62 


851569 


149.56419 


2.055025 


24.68 


3.50 


3.75 


-99 


3.75 


3 


no 


-20.84 


973691 


150.34480 


2.257383 


24.27 


5.08 


4.95 


-99 


4.95 


3 


no 


-21.55 


1200307 


150.28623 


2.424521 


23.34 


3.10 


1.18 


-99 


3.10 


3 


no 


-20.96 


1675730 


150.07668 


2.702570 


23.94 


3.10 


0.05 


-99 


3.10 


3 


no 


-21.07 


1234913 


150.11537 


2.363513 


23.06 


3.20 


3.46 


-99 


3.46 


3 


no 


-21.03 


1413401 


150.47971 


2.496383 


24.24 


3.12 


0.00 


-99 


3.12 


3 


no 


-20.65 


1897688 


150.07535 


2.828791 


23.78 


3.50 


0.61 


-99 


3.50 


3 


no 


-21.27 


502784 


150.40532 


1.910871 


24.09 


3.12 


2.62 


-99 


3.12 


^ 


no 


-20.44 


754632 


150.38487 


2.036528 


24.15 


3.12 


0.41 


-99 


3.12 


3 


no 


-20.79 


1464601 


150.21416 


2.475017 


23.19 


3.11 


0.81 


-99 


3.11 


3 


yes 


-21.71 


1498712 


149.76903 


2.573805 


23.41 


3.14 


0.24 


-99 


3.14 


3 


no 


-21.28 


1110682 


149.50595 


2.185332 


22.01 


3.28 


0.09 


-99 


3.28 


3 


no 


-22.71 


1781388 


149.50613 


2.642649 


24.94 


3.72 


1.40 


-99 


3.72 


3 


no 


-20.30 


265707 


150.70856 


1.698028 


23.86 


3.40 


0.17 


-99 


3.40 


3 


no 


-21.03 


593446 


150.05812 


1.794438 


24.44 


3.20 


3.16 


-99 


3.16 


3 


no 


-20.45 


644212 


149.53085 


1.953905 


24.28 


3.40 


0.21 


-99 


3.40 


3 


no 


-20.49 


1043168 


150.00752 


2.132265 


24.25 


3.40 


3.33 


-99 


3.33 


3 


no 


-20.83 


259947 


150.61942 


1.737010 


24.60 


3.10 


3.13 


-99 


3.13 


3 


no 


-20.46 


878986 


149.47784 


2.036321 


24.76 


4.10 


4.00 


-99 


4.00 


3 


no 


-20.61 


262808 


150.63403 


1.719672 


24.96 


4.60 


4.46 


-99 


4.46 


3 


no 


-20.84 


710344 


150.62828 


2.006204 


23.36 


3.60 


3.45 


-99 


3.45 


3 


no 


-22.06 


1384622 


150.59956 


2.534162 


23.29 


3.20 


0.53 


-99 


3.20 


3 


no 


-21.28 


702684 


150.56259 


2.060279 


24.34 


3.70 


3.70 


-99 


3.70 


3 


no 


-20.98 


1255908 


149.90977 


2.391198 


24.69 


3.90 


3.85 


-99 


3.85 


3 


no 


-20.85 


1723659 


149.77576 


2.704966 


17.94 


3.70 


3.77 


-99 


3.77 


2 


no 


-22.76 


428704 


149.55103 


1.764459 


24.59 


3.10 


3.20 


-99 


3.20 


2 


no 


-20.47 


462793 


150.73788 


1.884472 


24.97 


3.20 


1.05 


-99 


3.20 


2 


no 


-19.00 


1041149 


150.00462 


2.143163 


23.88 


3.22 


1.04 


-99 


3.22 


2 


no 


-22.82 


1257886 


149.97707 


2.378910 


24.71 


3.12 


0.21 


-99 


3.12 


2 


no 


-20.50 


1558001 


149.50140 


2.469177 


24.66 


5.12 


5.42 


-99 


5.42 


2 


no 


-22.67 


1309325 


149.56500 


2.362180 


24.65 


5.60 


5.28 


-99 


5.28 


2 


no 


-23.07 


1260818 


149.92712 


2.358468 


24.16 


3.20 


3.44 


-99 


3.44 


2 


no 


-21.06 


405715 


149.70297 


1.752171 


24.95 


3.10 


3.38 


-99 


3.38 


2 


no 


-20.27 


766066 


150.06465 


2.124108 


24.51 


5.00 


4.85 


-99 


4.85 


2 


no 


-21.93 


1631571 


150.51112 


2.665648 


24.61 


3.00 


3.82 


-99 


3.82 


2 


no 


-21.09 


1383442 


150.67668 


2.543255 


23.28 


3.20 


0.48 


-99 


3.20 


2 


no 


-21.09 


361093 


149.92401 


1.727810 


23.83 


3.33 


0.67 


-99 


3.33 


2 


no 


-20.86 


394142 


149.82733 


1.665374 


23.40 


3.12 


0.25 


-99 


3.12 


2 


no 


-21.53 


510917 


150.49745 


1.851680 


24.16 


4.60 


4.39 


-99 


4.39 


2 


no 


-21.82 


380027 


149.88055 


1.761200 


24.01 


3.60 


0.63 


-99 


3.60 


2 


no 


-21.02 


298002 


150.43706 


1.649305 


23.35 


3.30 


3.85 


3.89 


3.89 


2 


no 


-22.34 


203339 


149.51596 


1.609118 


23.76 


3.20 


3.49 


-99 


3.49 


2 


no 


-20.85 


964908 


150.42227 


2.150364 


23.97 


3.30 


0.28 


-99 


3.30 


2 


no 


-20.89 


1338373 


149.48837 


2.327380 


24.52 


5.10 


4.64 


-99 


4.64 


2 


no 


-21.66 


598191 


149.80116 


1.927920 


24.20 


4.60 


4.36 


-99 


4.36 


2 


no 


-21.75 


232585 


150.73338 


1.781481 


24.11 


3.20 


0.22 


-99 


3.20 


2 


no 


-20.91 


535642 


150.30261 


1.852066 


24.46 


3.40 


3.50 


3.84 


3.84 


2 


no 


-20.91 


1385272 


150.66348 


2.529267 


23.65 


3.30 


3.12 


-99 


3.12 


2 


no 


-21.16 


519634 


150.27715 


1.958373 


22.14 


3.40 


0.70 


-99 


3.40 


2 


no 


-22.61 


1110702 


149.55306 


2.185063 


23.20 


3.44 


3.46 


-99 


3.46 


2 


no 


-20.94 


1432719 


150.26303 


2.520890 


24.22 


3.30 


0.52 


-99 


3.30 


2 


no 


-20.69 


1508420 


149.83003 


2.507039 


24.78 


4.80 


0.63 


-99 


4.80 


2 


no 


-21.37 


184803 


149.60240 


1.580454 


23.42 


3.27 


0.41 


-99 


3.27 


2 


no 


-21.39 


316248 


150.25449 


1.694429 


24.49 


3.40 


0.42 


-99 


3.40 


2 


no 


-20.45 


1314059 


149.62923 


2.329145 


24.59 


3.40 


0.39 


-99 


3.40 


2 


no 


-20.61 


853930 


149.56989 


2.037617 


24.14 


3.12 


0.05 


-99 


3.12 


2 


no 


-20.55 
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TABLE 3 - Continued 



ID 


RA 
(J2000) 


DEC 
(J2000) 


Iab 
(AUTO) 


Zest 


Zqso 


Zspec 


z used 


Flag 


X-ray 


Mi4 50 


1615630 


150.55912 


2.769707 


24.95 


3.16 


3.35 


-99 


3.35 


2 


no 


-20.29 


1510107 


149.84480 


2.495914 


24.64 


3.70 


0.64 


-99 


3.70 


2 


no 


-20.70 


840968 


149.81004 


2.124751 


23.65 


3.80 


0.63 


-99 


3.80 


2 


no 


-21.93 


418417 


149.56200 


1.662382 


24.05 


5.30 


0.96 


-99 


5.30 


2 


no 


-22.99 


1596481 


150.67618 


2.732485 


24.85 


3.40 


3.12 


-99 


3.12 


2 


no 


-19.76 


1045616 


149.79565 


2.283589 


24.80 


4.10 


4.21 


-99 


4.21 


2 


no 


-20.93 


1385228 


150.64104 


2.529159 


23.51 


3.30 


0.62 


-99 


3.30 


2 


no 


-21.26 


1552773 


149.47151 


2.514272 


23.66 


4.14 


0.21 


-99 


4.14 


2 


no 


-21.19 


1893339 


150.17361 


2.856197 


24.49 


4.26 


4.22 


-99 


4.22 


2 


no 


-20.97 


697377 


150.59901 


2.099155 


23.10 


3.22 


0.29 


-99 


3.22 


2 


no 


-21.47 


1943641 


149.78291 


2.824640 


24.26 


5.00 


5.02 


-99 


5.02 


2 


no 


-22.07 


1739756 


149.60091 


2.764447 


24.34 


5.80 


5.45 


-99 


5.45 


2 


no 


-23.48 


851974 


149.58417 


2.046672 


24.32 


3.40 


3.66 


-99 


3.66 


2 


no 


-20.98 


1665352 


150.12608 


2.767497 


24.80 


6.00 


2.37 


-99 


6.00 


2 


no 


-22.25 


1401595 


150.46024 


2.577492 


24.55 


3.28 


0.21 


-99 


3.28 


2 


no 


-20.43 


784661 


150.10434 


2.001693 


24.21 


3.30 


0.27 


-99 


3.30 


2 


no 


-20.82 



